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1.  INTRODUCTION 


The  growing  emphasis  on  energy  conservation  has  led  to  accelerated  research 
on  the  development  of  efficient  equipment  for  heat  exchange  processes.  Today, 
a  large  share  of  this  research  is  directed  toward  techniques  for  augmenting  heat 
transfer  [1.1-1.3].  The  ultimate  purpose  of  such  techniques  may  be  to  improve 
the  thermodynamic  performance  of  future  or  existing  heat  exchangers,  to  reduce 
their  size,  thereby  reducing  capital  and  operating  costs,  and  possibly,  to  prevent 
excessive  temperatures  and  failure  in  systems  where  the  heat  generation  rates 
are  difficult  to  tolerate. 

The  task  of  evaluating  the  merit  of  a  proposed  augmentation  technique  may 
be  as  important  as  developing  the  technique  itself.  Experts  agree  that  evaluating 
the  impact  of  these  techniques  is  a  difficult  problem  fl.1-1.3),  partly  due  to 
their  great  diversity  and  partly  due  to  the  absence  of  a  common  conceptual  basis 
for  evaluation.  The  present  study  was  motivated  by  the  real  need  for  a  univer¬ 
sal,  thermodynamically  solid,  basis  for  comparing  the  worth  of  all  augmentation 
techniques. 

In  general,  the  objective  of  all  heat  transfer  augmentation  techniques  is 
to  enhance  the  thermal  performance  by  increasing  the  surface  heat  transfer 
coefficient  relative  to  the  heat  transfer  coefficient  which  characterizes  the 
standard  (untouched)  surface.  At  the  same  time,  the  increase  in  heat  transfer 
coefficient  should  be  accompanied  by  only  a  small  increase  in  the  fluid  pumping 
power  loss  needed  for  the  heat  exchange  job.  These  objectives  underlie  the 
fundamental  trade-off  in  the  development  of  any  heat  exchange  equipment,  namely, 
a  design  modification  which  improves  the  thermal  contact  will  most  likely  cause 
a  parallel  increase  in  the  mechanical  power  dissipated  in  the  apparatus,  for 
example  the  power  needed  for  pumping  the  working  fluid  through  the  passages  of 
a  heat  exchanger.  Given  this  trade-off,  it  is  important  to  know  in  advance 
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which  augmentation  technique  will  realistically  lead  to  an  improved  overall 
performance.  In  addition  to  finding  out  whether  an  improvement  is  feasible, 
it  is  important  to  know,  in  quantitative  thermodynamic  terms,  the  size  of 
the  improvement. 

This  report  summarizes  the  main  results  of  a  two-year  study  of  the  impact 
of  heat  transfer  augmentation  on  the  thermodynamic  performance  of  heat  exchange 
equipment.  The  study  relied  on  the  concept  of  thermodynamic  irreversibility 
and  entropy  generation  as  fundamental  criterion  for  evaluating  and,  eventually, 
minimizing  the  waste  of  usable  energy  (exergy)  in  energy  systems.  The  extreme 
importance  of  this  thermodynamic  concept  in  the  world  of  energy  engineering 
was  summarized  most  recently  at  the  1979  Second  Law  of  Thermodynamics  Workshop 
[l.4].  It  has  also  been  shown  [l. 3, 1.5, 1.6]  that  the  ultimate  function  of  all 
augmentation  techniques  is  to  inhibit  the  production  of  entropy,  thus  paving 
the  road  toward  thermodynamically  efficient  equipment  for  heat  exchange. 

The  work  summarized  in  this  report  is  organized  into  three  distinct  parts. 
The  first  part  (section  2)  contains  a  fundamental  investigation  of  the  entropy 
generation  mechanism  in  elementary  heat  transfer  configurations.  This  tutorial 
section  establishes  the  foundation  for  the  second  law  analysis  of  actual  aug¬ 
mentation  techniques,  reported  in  sections  3  and  4.  Section  3  describes  in 
quantitative  terms  the  irreversibility  minimization  potential  of  some  of  the 
most  common  techniques,  namely,  swirl-flow  devices  and  roughened  surfaces. 
Section  4  considers  the  entropy  generation  penalty  associated  with  the  use  of 
extended  surfaces  (fins)  in  convective  heat  transfer.  This  section  shows 
also  how  the  fin  geometry  can  be  chosen  such  that  the  fin  performs  its  heat 
transfer  duty  with  minimum  generation  of  entropy. 


Due  to  space  limitations,  the  present  report  reviews  only  a  segment  of 
the  research  undertaken  during  this  two-year  study.  For  a  complete  and  detailed 


account  of  all  the  topics  covered  by  this  study,  the  reader  is  encouraged  to 
consult  References  [1.7-1.12]. 
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2.  ENTROPY  GENERATION  IN  CONVECTIVE  HEAT  TRANSFER 

2 . I  Local  Rate  of  Entropy  Generation 

Consider  the  two-dimensional  infinitesimal  fluid  e.  ement  dxdy  shown 
schematically  in  Fig.  2.1.  The  fluid  element  is  part  of  a  considerably  more 
complex  convective  heat  transfer  picture.  However,  for  the  scope  of  this 
study,  wo  regard  the  element  as  an  open  thermodynamic  system  subjected  to  mass 
fluxes,  energy  transfer  and  entropy  transfer  interactions  through  a  fixed  con¬ 
trol  surface.  The  element  size  is  small  enough  so  that  the  thermodynamic  state 
of  the  fluid  inside  the  element  may  be  regarded  as  uniform  (independent  of 
position) .  However,  the  thermodynamic  state  of  the  small  fluid  element  may 
change  with  time. 

For  this  study,  we  limit  our  attention  to  incompressible  fluids  without 
internal  heat  generation.  In  such  cases,  the  expression  for  the  volumetric  rate 
of  entropy  generation  reduces  to  (2.1] : 


As  Wf>  might  have  expected,  the  irreversibility  indicator  S'"  contains  two  addi¬ 
tive  parts,  one  due  to  conduction  in  the  presence  on  non-zero  temperature 
gradients,  the  other  accounting  for  viscous  dissipation  of  mechanical  power 
throughout  the  flow.  In  the  second  term  of  equation  (2.1),  the  factor  in 
brackets  represents  the  viscous  dissipation  function  $  for  two-dimensional 
compressible  flow. 

The  local  entropy  generation  rate  formula  (2.1)  can  be  derived  in  straight¬ 
forward  fashion  by  performing  an  entropy  transfer  accounting  around  the 
infinitesimally  small  element  illustrated  in  Fig.  2.1.  The  entropy  transfer 


Entropy  generation  analysis  for  an  infinitesimal  element  dxdy  in 
convective  heat  transfer. 
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Left:  entropy  generation  profiles  N  ,  inside  a  smooth  tube  with 

constant  wall  heat  flux  in  laminar  flow.  Right:  optimum  Reynolds 
number  and  corresponding  minimum  entropy  generation  number  for 
forced  convection  heat  transfer  in  a  round  tube. 
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t.o  and  from  tlio  (dxdy)  system  is  associated  witii  heat  transfer ,  q  and 

x 

q  ,  as  well  as  with  mass  Fluxes,  pv  and  pv  .  In  the  interest  of  brevity 
y  x  y  1 

we  do  not  show  this  derivation,  urginq  the  reader  to  consult  any  irreversible 
thermodynamics  book  for  additional  details.  Alternate  versions  of  formula 
(2.1)  corresponding  to  coordinate  systems  other  than  the  two-dimensional  carte¬ 
sian  system  of  Pig.  2.1  may  be  found  in  [2.2]. 

In  many  heat  transfer  problems  it  is  often  possible  and  convenient  to 
neglect  the  viscous  dissipation  term  n<J>  in  the  equation  for  energy  conservation. 


/  38  30  \ 

1  o2e  a2e  \ 

\  x  Ox  +  Vy  3y  ) 

=  k 

'  0x2  +  3y2  / 

This  is  particularly  the  case  in  heat  transfer  through  gases  at  subsonic  velo¬ 
cities.  Tiie  dimensionless  group  which  expresses  the  magnitude  of  the  dissipation 
energy  term  relative  to  the  conduction  energy  term  in  (2.2)  is  [2.3] 


be  Pr 


u*2  v. 

c  0*  a 
P 


(2.3) 


where  u*  and  G*  are  the  characteristic  fluid  velocity  and  temperature 
difference  for  the  convective  heat  transfer  problem.  Thus,  in  many  engineering 
problems,  we  find  Ec  Pr  <<  1  . 

Consider  now  the  same  question  relative  to  expression  (2.1).  Under  what 

conditions  is  the  viscous  dissipation  contribution  to  S'"  negligible?  If  we 

regard  expression  (2.1)  as  the  sum  S'"  =  S'"  ,  +  S'"  .  ,  then,  in 

conductive  viscous 

an  order  of  magnitude  sense. 


viscous 


Ec  Pr 


conductive  > 


(2.4) 
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llere,  i  '  ' 1  * /T *  ,  where  T*  is  the  absolute  temperature  characteristic  to 

the  problem  at  hand.  The  dimensionless  tempera! ure  difference ,  t  ,  is  always 
an  important  dimensionless  parameter  in  second  law  analyses  of  heat  transfer 
problems.  With  the  exception  of  applications  at  cryogenic  temperatures,  the 
temperature  difference  number  c  is  generally  much  smaller  than  unity, 

•i  <<-  1  . 

It  is  now  clear  that  the  energy  argument  by  which  viscous  dissipation  is 
neglected  in  (2.2)  ,  i.e.,  Ec  Pr  <<  1  ,  has  no  bearing  on  the  question  of 
negligible  viscous  contribution  to  the  local  rate  of  irreversibility  production. 
It.  is  then  possible  to  encounter  situations  where,  although  the  energy  equation 
can  be  simplified  according  to  Ec  Pr  <<  1  ,  S'"  is  in  fact  dominated  by 

viscous  effects.  This  is  the  limit  in  which  t  is  very  small,  small  enough 
so  that  Ec  Pr/x  >  1  . 

Below  we  examine  a  series  of  important  convective  heat  transfer  configura¬ 
tions  in  an  effort  to  illustrate  the  coupling  of  viscous  and  conductive  effects 
in  the  makeup  of  S'"  .  In  the  process  wc  will  study  the  spatial  distribution 
of  irreversibility,  pointing  out  those  flow  features  which  act  as  concentrators 
(sources)  of  entropy  generation  S’"  . 


2 . 2  Forced  Convection  in  a  Round  Tube 

Laminar  Flow.  Consider  the  Poiseuille  flow  through  a  round  tube  with  uni¬ 
form  heat  flux  q"  around  its  circumference  (see  insert  in  the  left  side  of 
Fig.  2.2).  The  velocity  and  temperature  profiles  for  this  flow  are  particularly 
simple  [2.2] : 
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(2.0) 
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x  r  v 
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Ttie  equation  for  S'"  in  the  cylindrical  geometry  of  Fig.  2.2  is 


S'  = 


which  yields 


S'"  = 


q"2 

kT? 


(2R  -  R3)  2  + 


Pe2  J 


4pv' 


Tr 


max 

1“ 


R  = 


(2.9) 


(2.10) 


liquation  (2.10)  is  the  entropy  generation  profile  in  the  pipe  cross-section. 
Together  with  the  velocity  and  temperature  profiles,  the  entropy  generation 
profile  completes  the  thermodynamic  description  of  the  convective  heat  transfer 
phenomenon. 

It  is  convenient  to  nondimensionalize  expression  (2.10)  and  define  the 
local  entropy  generation  number  N 


kT 


S’" 


=  S" 
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j"  2 
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R3)2 
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Pe2 


J±£-l£  r2 


(2.11) 


Hi  re  we  made  the  assumption  that  the  temperature  variation  over  the  pipe 

cross-section  is  negligible  compared  with  the  absolute  temperature,  hence 

T  -  T  where  T  is  a  characteristic  (reference)  absolute  temperature.  The 
o  o 

irreversibility  figure  Ns„,  depends  on  radial  position,  on  Pe  and  on  the 
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group  Ec  Pr/x  which,  as  discussed  in  the  preceding  section,  determines  the 
relative  importance  of  viscous  effects.  The  Peclet  number  Pe  governs  the 
importance  of  irreversibility  associated  with  conduction  in  the  axial  direction. 
We  see  that  for  Pe  <  4  the  axial  conduction  contribution  dominates  the  radial 
conduction  effect. 

Tile  left  side  of  Fig.  2.2  displays  a  family  of  entropy  generation  profiles 

in  the  pipe  cross-section,  for  cases  where  the  axial  conduction  effect  is 

negligible  (Pe  >>  4) .  The  value  of  Ec  Pr/x  increases  gradually  to  the  point 

where  viscous  effects  dominate  N  m  .  In  all  cases,  the  pipe  wall  region  acts  as 

a  strong  concentrator  of  irreversibility.  When  Ec  Pr/x  =  0  ,  the  maximum  N  ,,, 

1/2 

occurs  inside  the  fluid  at  R  =  (2/3)  ,  the  same  place  where  due  to  the  wall 

curvature  effect  the  maximum  radial  temperature  gradient  is  located.  As  viscous 
effects  take  over,  the  point  of  highest  Ng,,,  migrates  toward  the  wall  and,  for 
Ec  Pr/x  >  1/4  ,  it  coincides  with  the  wall. 

The  rate  of  entropy  generation  over  the  entire  tube  cross-section  is 
obtained  by  integrating  expression  (2.10), 


S ' 


2n  r 

o 


R  d  R  . 


(2.12) 


Neglecting  axial  conduction,  Pe  >>  4  ,  result  (2.12)  can  be  written  as 


11  g' 2  8_  ym2 

48tt  kT2  it  Tp 2  r  **  ' 


(2.13) 


Hero  m  is  the  mass  flow  rate  through  the  tube,  while  q’  is  the  heat  transfer 
rate  per  unit  length,  2-iTr^q"  .  Once  more,  the  irreversibility  production  S' 
appears  as  the  sum  of  two  effects,  heat  transfer  in  the  direction  of  a  finite 
temperature  gradient  and  fluid  friction.  It  should  be  noted  that  the  heat  transfer 
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contribution  to  S'  ,  the  first  term  in  equation  (2.13),  is  fixed  as  soon  as 
the  heat  transfer  rate  per  unit  length  q'  (the  heat  transfer  duty  of  the 
tube)  is  specified.  We  return  to  this  observation  later  in  this  section  when 
we  address  the  question  of  optimum  tube  radius  for  minimum  irreversibility 
and  fixed  q'  . 

Turbulent  Flow.  Exact  analytical  solutions  for  the  turbulent  velocity 
and  temperature  fields  in  the  tube  cross-section  are  not  available.  Therefore, 
one  cannot  evaluate  the  rate  of  entropy  generation  at  any  point  in  the  pipe 
cross-section,  as  done  through  equation  (2.10)  for  laminar  flow.  However,  one 
can  still  evaluate  the  rate  of  entropy  generation  averaged  over  the  tube  cross- 
section  by  using  the  integral  result  developed  by  Bejan  [2.4]  for  irreversibility 
in  a  duct  of  arbitrary  geometry. 

Using  the  present  notation  for  forced  convection  through  a  tube,  the 
integral  result  is  [2.4] 


1  q'2 

TT  Nu  kT2 


+ 


m3f 

p?-T  r  5 
o 


(2.14) 


which  shows  how  S'  can  be  evaluated  based  on  average  heat  transfer  (Nu)  and 
fluid  friction  (f)  information.  As  one  might  expect,  the  laminar  flow  expres¬ 
sion  (2.13)  is  only  a  special  case  of  the  more  general  result  (2.14)  since, 

in  laminar  flow,  Nu  =  48/11  and  f  =  8  tr  y  r  /m  .  Unlike  in  laminar  flow, 

o 

the  heat  transfer  contribution  to  S '  is  not  necessarily  constant  when  the 
heat  transfer  duty  q'  is  specified. 

Optimum  Tube  Radius  for  Minimum  Irreversibility.  In  an  application  in 
which  the  heat  transfer  duty  q'  and  the  mass  flow  m  are  already  specified, 
it  is  possible  to  select  an  optimum  tube  radius  which  insures  the  minimum  rate 
of  entropy  generation  in  the  heat  transfer  device.  This  design  optimization 


¥ 


1 
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procedure  is  described  best  by  placing  the  irreversibility  rate  expression 

(2.14)  in  dimensionless  form.  We  define  the  entropy  generation  number  ' 

as  the  ratio  S' /S'  where  S'  is  tin;  actual  entropy  generation 

AT, laminar 

rate  qiven  by  (2.14)  and  S'  ,  is  the  first  term  appearing  in  (2.13). 

AT, laminar 

We  commented  earlier  that  S!  ,  .  is  constant  when  the  heat  transfer 

At, laminar 

rate  q'  is  specified,  lienee.  S'  ,  assumes  the  role  of  characteristic 

1  1  AT, laminar  - 

rate  of  entropy  generation  in  the  system  of  interest. 

The  entropy  generation  number  is  therefore 


N 


S' 


48 

11 


Nu 


-1 


t' 


Re‘J 


(2.11) 


showing  that  the  duct  irreversibility  depends  primarily  on  two  dimensionless 
groups,  the  Reynolds  number  based  on  tube  diameter  Re  and  the  "duty" 
parameter  , 


Re 


2m 

PUro  ; 


B 

o 


q'm 


h5/2(kT)1/2 


(2.16,17) 


With  q'  ,  m  and  working  fluid  specified,  the  task  of  finding  the  optimum  radius 
for  minimum  S"  is  equivalent  to  minimizing  the  N<. '  expression  (2.15)  with 
respect  to  Re  ,  subject  to  a  specified  constant  .  This  procedure  is 
straightforward,  therefore  only  the  final  results  are  given  here. 

(•’or  laminar  flow,  using  Nu  =  48/11  and  f  =  16/Re  in  expression  (2.15), 
the  minimization  procedure  yields 

Re  _  =  0  ;  N  , mln  =  1  (2. 18, Id) 

opt  S 

In  engineering  terms,  this  result  implies  that  the  selected  tube  radius  r 
must  be  large  enough  so  that  the  rate  of  entropy  generation  is  strongly 


dominated  by  the  contribution  due  to  heat  transfer  across  a  finite  temperature 
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dif lor^ni'i  .  In  other  words,  based  on  expression  (2.15), 


.’4tt 

11 


Re‘*B 


(2.20) 


For  turbulent  flow,  the  N^. ,  expression  (2.1r>)  has  a 
Substituting  Nu  -  0.023  Pr^’^Re^"^  and  f  =  0.046  Re 
dif  l'erentiat  mg  with  respect  to  Re  yields 


unique  minimum, 
into  (2.15)  and 


Re 

opt 


-0.071  -0.356 

023  Pr  B 

o 


(2.21) 


and 


N  .  . m  L "  ^  1- 


-0.343  -0.2o6 

Pr  B 

o 


(2.22) 


Expressions  (2.21,  22)  nave  been  summarized  in  the  right  hand  side  of  Fig.  2.2 
for  two  discrete  values  of  Prandtl  number.  As  the  aggregate  duty  parameter 


B  increases,  we  see  that  the  optimum  tube  radius  decreases  (Re  increases) 
■■  >  opt 

and  thi.-  minimum  entropy  generation  number  N  mln  decreases  also. 


2.  >  boundary  Layer  Over  Flat  Plate1 

bamrnar  Flow  over  Isothermal  Plate.  Consider  now  the  development  of 
laminar  momentum  and  thermal  boundary  layers  along  a  flat  plate.  The  situation 
is  shown  schematical ly  in  the  horizontal  plane  of  the  isometric  drawing  of 
Fig.  2.3.  At  some  distance  from  the  solid  wail  the  fluid  velocity  and  temperature 
are  uniform,  and  T  .  The  wall  temperature  is  constant,  T  . 

x,oo  or.  0 

The  study  of  the  velocity  and  temperature  fields  in  the  vicinity  of  the 
plate  constituted  the  subject  of  numerous  investigations  [2.5].  The  purpose 
of  this  section  is  to  examine  the  distribution  of  entropy  generation  in  the 


boundary  layer.  For  this  we  rely  on  solutions  available  in  the  literature  for 
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v  (x,y)  and  0(x,y)  in  laminar  flow, 
x 

The  task  of  evaluating  the  entropy  generation  profile  S'"  is  simplified 
greatly  if  we  restrict  the  discussion  to  the  case  Pr  =  1  for  which  the 
Blasius-Pohlhausen  solution  [2.6]  reduces  to 


0  Vx  df 


0  v  dg 

oo  X,  <"> 


(2.23,24) 


The  similarity  variable  g  equals  y[v^  M/ (vx) 3  <  while  f(g)  is  the 

function  tabulated  by  Howarth  [2.7].  Neglecting  the  irreversibility  terms 
associated  with  velocity  and  temperature  gradients  in  x  direction,  we  find 


S 


•  II 


kO  7v 
°°  x,1" 

T  ?vx 
o 


uv 


(f") 


X,* 


T  VX 
o 


(f")2 


(2.25) 


The  local  entropy  generation  number  is 


Ec  Pr  \ (f")2 

+  -  1— - ,  Pr 

r  /Re 


1  , 


(2.26) 


wlii  re  Re  is  defined  as  v  x/v  . 

X  X,” 

The  complex  dependence  of  S'"  on  both  x  and  y  is  shown  in  Fig.  2.3. 

The  three-dimensional  display  was  done  in  terms  of  (x  v  /v)  and 

x 

(y  v  /v)  in  the  horizontal  plane,  and  N_tll/(1  +  EcPr/t)  in  the  vertical 
x,  ^  S 

direction.  It  is  evident  that  the  irreversibility  effects  are  limited  to  the 
boundary  layer.  Regarding  the  y  dependence  of  N  ,  the  entropy  generation 
rate  is  highest  reaching  a  peak  at  the  solid  wall.  The  longitudinal  variation 
of  S'"  is  as  1/x  ,  indicating  that  like  all  gradients  in  this  boundary  layer 
solution  S'"  blows  up  at  the  origin.  The  viscous  effect  again  scales  up  as 


Ec  Pr/t  . 


flow  direction 


F.ntropy  generation  surface  for  laminar  boundary  layer  flow  and 
heat  transfer  along  a  flat  plate. 


Optimum  Reynolds  number  and  corresponding  minimum  entropy  generation 
number.  Left:  boundary  layer  flow  over  a  flat  plate.  Riqht: 
single  cylinder  in  gaseous  cross-flow. 
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Integrating  (2.25)  across  the  boundary  layer  we  can  calculate  the  rate 
of  entropy  generation  per  unit  area  of  flat  plate, 


kv  6 


v  00  cx>  /  EC  Pr  \ 

S"  =  S-  dy  =  0.25  -  ( 1  +  ) 


(2.27) 


with  the  corresponding  entropy  generation  number  defined  as 


v  =  s"  k7~r? 

V  oo  00 


(2.28) 


Finally,  by  analogy  with  the  Nusselt  number  nomenclature  for  boundary 
layer  heat  transfer,  we  integrate  (2.27)  in  the  x  direction  to  find  the 
total  rate  of  entropy  generation  produced  by  boundary  layer  flow  and  heat 
transfer  ov<  r  a  length  L 


S"  dx  =  0.50 


I ,  Ec  Pr  \ _ 

(i  +  “T- )nei 


(2.29) 


wtr  re  K<-[  -  v^  _E/v  .  The  overall  entropy  generation  number  based  on  S  ’  is 


-  O.b0  (  i  ♦  ES-EE. ), 


(2.30) 


It  is  wortn  mentioning  that  in  general  tlie  Prandtl  number  will  have  an 
additional  effect  on  the  relative  importance  of  viscous  and  conductive  effects 
in  t.h<-  constitution  of  S’"  ,  S"  and  S’  .  It  is  easy  to  show  that  when  Pr  ^  1 
in.  viscous  effects  scaie  as 


Ec  Pr  |  jr j' .  Ec  Pr1/j 


(2.31) 
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where  6  and  5^  are  the  velocity  and  temperature  boundary  layer  thicknesses. 

In  writing  (2.31)  we  made  use  of  the  approximation  6/6,^  *  Pr^"*  .  The 
significance  of  (2.31)  is  that  when  the  two  thicknesses  are  not  equal,  the 
thinner  layer  exhibits  larger  gradients  thereby  enhancing  its  contribution  to 
the  total  S'"  figure.  According  to  (2.31),  viscous  effects  are  more  likely  to 
play  a  role  in  low  Prandtl  number  fluids  where  6  <<  6^  . 

Laminar  and  Turbulent  Flow  over  Constant  Heat  Flux  Plate.  To  study  the 
generation  of  entropy  in  a  turbulent  boundary  layer  one  has  to  rely  on  an  integral 
method  which  takes  into  account  the  heat  transfer  and  fluid  friction  characteristics 
of  the  flow  in  an  overall  manner,  as  Nusselt  number  and  friction  factor  informa¬ 
tion.  Consider,  for  example,  a  flat  plate  of  negligible  thickness  suspended  in 

a  uniform  flow  field,  parallel  to  the  flow  velocity  v  .  The  heat  flux  q" 

x,°° 

over  the  plate  surface  is  uniform.  Consider  also  a  control  surface  which  sur¬ 
rounds  the  plate  of  finite  length  L  at  a  large  enough  distance  through  portions 
of  external  flow  in  which  the  fluid  motion  is  nearly  uniform  and  the  temperature 
nearly  constant,  T^  .  Regardless  of  whether  the  boundary  layer  is  turbulent  or 
laminar,  the  entropy  generation  rate  in  one  half  of  the  control  volume  (i.e., 
for  one  side  of  the  plate)  is  given  by 


S' 


2  T 

03 


L 


0 


(2.32) 


where  h  and 


x 

coefficient, 


C 


f  ,x 


are  the  local  heat  transfer  coefficient  and  skin  friction 


h 

x 


q" 

T  (x)  -  T 
o  °° 


C 


f  ,x 


2  T 

o 


pv 

X,” 


(2.33,34) 
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In  the  above  definitions  T  (x)  is  the  wall  temperature,  r  the  wall  shear 

o  o 

stress  and  q"  the  uniform  heat  flux.  Note  that  unlike  equation  (2.29)  in  the 
preceding  sub-section  where  we  considered  the  laminar  boundary  layer  over  an 
isothermal  flat  plate,  the  S'  expression  (2.32)  refers  to  a  uniform  heat 
flux  situation. 

Expression  (2.32)  is  the  result  of  an  entropy  flux  accounting  around  the 

control  volume,  analysis  omitted  here  due  to  space  limitations.  Like  all 

entropy  generation  results  for  forced  convection  heat  transfer.  S'  consists 

of  two  additive  parts,  one  due  to  heat  transfer  across  the  [T  (x)  -  T  ] 

o  00 

temperature  difference,  the  other  being  associated  with  the  total  friction  drag 
force  exerted  by  the  fluid  on  the  plate.  Below,  we  use  result  (2.32)  to  deter¬ 
mine  the  optimum  plate  length  L  which  yields  the  minimum  rate  of  entropy 
generation  in  a  heat  transfer  application  in  which  the  uniform  flow  velocity 
v^  ^  and  the  total  heat  flux  qr  =  j  ^q"  dx  are  specified.  We  do  this  by 
first  substituting  appropriate  correlations  for  h  and  C  into  expression 

X  I  »x 

(2.32)  and  solving  the  equation  9S'/3L  =  0  . 

For  laminar  flow,  the  local  skin  friction  coefficient  is  =  0.664  Re 

f  ,x  x 

while  the  local  heat  transfer  coefficient  is  given  by  h  x/h  =  0.332  Pr^3  Re 

x  x 

[2.8] .  Writing  q'  =  q"L  for  the  total  heat  transfer  rate  from  plate  to  fluid 
over  the  plate  length  L  ,  the  entropy  generation  number  Ng ,  becomes 


-1/2 


S  ’ 


kT  z 

oo  —1/3  —1/2  1/2  —2 

=  S'  — TJ-  =  2.008  Pr  '  Re„  '  +  0.664  Re  /  B 

n  '  ^  T.  T. 


(2.35) 


Here,  Re  and  B  are  the  Reynolds  number  based  on  L  and  the  "duty" 
h 

parameter,  respectively. 


Re  = 

L 


v  L 

X,” 


B  = 


v  (vikT) 

X,“ 


1/2 


(2.36,37) 
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Tliu  optimum  plate  length 
tion  at  constant  q'  and 


Rer 

L,opt 
v  is 


yielding  the  minimum  rate  of  entropy  genera- 


x,1” 


Re  =  3.024  Pr 

L,opt 


-1/3 


=  2.309  Pr"1/6  B'1 


(2.38,39) 


For  turbulent  flow  we  use  a  similar  set  of  correlations  for  friction  and 
heat  transfer  =  0.0576  Re  ^ ^  and  h  x/k  =  0.0296  Pr3^3  Re  0,0  [15]. 

f,X  XX  X 

We  also  assume  that  the  laminar  layer  which  precedes  the  turbulent  boundary 
layer  is  much  shorter  than  the  plate  length  L  .  Substituting  these  correla¬ 
tions  into  the  entropy  generation  result  (2.32)  yields 


kT  ^ 

«  -1/3  8/5  -8/5  -2 

NS'  =  S'  ~q77'  =  28'15  Pr  ReL  +  0,036  ReL  B  • 


Re  ^  =  64.31  Pr“5/24  B574,  =  2.013  Pr'1/6  b’1  . 

L,opt  S" 


(2.40) 


(2.41,42) 


The  optimum  plate  length  and  the  resulting  minimum  entropy  generation 
number  prescribed  by  these  equations  are  shown  on  the  left  side  of  Fig.  2.4. 

The  discontinuity  illustrated  by  dashed  lines  corresponds  to  the  transition 
region,  cases  in  which  the  laminar  and  turbulent  portions  of  the  boundary 
layer  are  of  comparable  lengths.  The  trends  are  similar  to  those  presented  in 
the  right  graph  of  Fig.  2.2:  the  higher  the  duty  parameter  B  ,  the  higher  the 
optimum  Reynolds  number  (plate  length)  and  the  lower  the  minimum  entropy 
generation  number.  In  other  words,  if  the  total  heat  transfer  rate  q1  is 
constant  and  the  flow  velocity  increases,  the  optimum  plate  length  decreases 
and  the  minimum  attainable  entropy  generation  rate  increases. 

These  conclusions  are  applicable  to  the  optimization  of  local  flow  geometry 
in  plate-fin  surfaces  for  compact  heat  exchangers,  as  shown  in  Chapter  4  of  this 


report . 
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2 . 4  Single  Cylinder  in  Cross-Flow 

The  heat  transfer  between  a  cylindrical  surface  and  a  fluid  flow  normal 
to  tlu'  cylinder  axis  is  one  of  the  most  frequent  heat  transfer  configurations 
encountered  in  actual  engineering  equipment  for  heat  exchange  processes.  It 
is  appropriate  to  examine  here  the  thermodynamic  irreversibility  introduced 
by  this  configuration.  Due  to  inherent  similarities  with  the  plate  of  finite 
length  in  parallel  boundary  layer  flow,  example  concluded  in  the  preceding 
section,  the  case  o'"  a  single  cylinder  in  cross-flow  will  be  summarized  very 
briefly  and  only  the  key  analytical  results  will  be  given. 

Consider  a  cylinder  of  diameter  D  with  uniform  surface  heat  flux  q" 
in  a  cross-flow  of  uniform  velocity  v^_  ^  and  temperature  .  An  entropy 
generation  analysis  of  the  flow  region  affected  by  heat  and  momentum  transfer 
from  the  cylinder  yields  the  total  entropy  generation  rate 

"  v  °oFn 

S'  =  “Sr  4  [TW  -  T  ]  dZ  +  ~X~° .  (2.43) 

T  2  T  «>  T 

OO  J  CO 


In  this  expression,  !  is  the  curvilinear  coordinate  around  the  cylindrical 
surface  and  T(£)  is  the  local  surface  temperature.  It  is  again  assumed  that 
tin.*  temperature  differences  are  much  smaller  than  the  characteristic  absolute 
temperature  of  the  medium  (T(£)  -  T  ]  <<  T  .  F  is  the  force  (drag)  per 

CO  Q 

unit  length  exerted  by  flow  on  cylinder,  force  calculated  from  drag  coefficient 
experimental  information,  CQ  , 


F  /D 
D 

:  “  pv  2/2 
x,°° 


(2.44) 


Replacing  the  lino  integral  in  (2.43)  with  the  average  wall-fluid  temperature 


difference  times  nu  ,  and  expressing  the  average  temperature  difference  in  terms 
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of  tilt:  average  Nusselt  number,  T(Jt)  -  T  •-=  q"D/ (k  Nu)  ,  yields, 
J  av  «■ 


S’ 


q‘2  _1_  +  1 

nkT  2  Nu  2D 

03 


Re 


T 

oo 


(2.45) 


Here,  we  used  the  drag  coefficient  to  replace  F^  in  equation  (2.43).  We 
also  wrote  q‘  =  TiDq"  for  the  total  heat  transfer  rate  per  unit  length  of 
cylinders.  Finally,  expression  (2.45)  is  put  in  dimensionless  form  defining 
the  entropy  generation  number 


kT 


Ns,  =  S' 


—  +  T  Cn  Re  B-2 
irNu  2  D 


(2.46) 


where  B  is  a  duty  parameter  which  has  the  same  form  as  in  equation  (2.37). 

As  in  the  heat  transfer  configurations  examined  earlier,  we  can  use  the 

entropy  generation  number  formula  (2.46)  to  determine  which  flow  geometry 

(cylinder  diameter  D)  is  best  for  minimizing  the  thermodynamic  losses  associated 

with  th"  heat  transfer  process.  For  an  application  in  which  q',  v  and 

x,°° 

the  fluid  are  known,  the  optimization  procedure  amounts  to  minimizing  expression 

(2.46)  with  respect  to  Re  ,  subject  to  constant  B  .  The  results  of  this 

optimization  procedure  are  shown  as  Re  and  Nmtn  on  the  right  side  of 

opt  S ' 

Fig.  2.4,  a  plot  qualitatively  similar  to  the  left  graph  obtained  for  a  flat 
plate.  Tlie  right  graph  of  Fig.  2.4  was  constructed  based  on  equation  (2.46) 
coupled  with  Hilpert's  average  Nusselt  number  correlation  for  gas  flow  [2.9] 
and  with  Kisner's  presentation  of  drag  coefficient  [2.10]. 


2 . 5  Laminar  Forced  Convection  in  the  Entrance  Region  of  a  Flat  Rectangular  Duct 
A  frequent  flow  configuration  which  embodies  the  characteristics  of  the 
tirst  two  flow  examples  discussed  here  is  sketched  in  the  insert  of  Fig.  2.5. 

In  th<  entrance  region  of  a  parallel-plate  duct,  temperature  and  velocity  boundary 
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layers  dovriop  simultaneously  along  both  walls,  gradually  filling  the  duct  and 
leading  to  the  well-known  fully-developed  laminar  regime.  There  have  no.  n  many 
studies  of  the  heat  transfer  and  velocity  problems  associated  with  this  basic 
configuration,  as  indicated  recently  by  Bhatti  and  Savory  [2.11J.  In  this 
section  wo  address  the  irreversibility  problem,  once  again  relying  on  published 
olutions  for  the  temperature  and  velocity  field  in  the  duct. 

Sparrow  [2.12]  analyzed  the  boundary  layer  development  in  the  entrance 
region  using  the  Karman-Pohlhausen  integral  technique.  The  Sparrow  solution 
is  unique,  considering  its  simplicity  versus  the  good  agreement  between  its  pre¬ 
dictions  and  experimental  findings.  In  what  follows  we  use  this  solution 
and  the  S'"  formula  (2.1)  to  illustrate  the  distribution  of  irreversibility 
in  the  channel. 

The  ultimate  goal  of  the  irreversibility  analysis  is  to  discover  how  the 
thermodynamic  losses  due  to  entrance  effects  compare  with  similar  losses  in  the 
lully-developed  section  of  the  duct.  In  other  words,  what  length  of  fully- 
devolopod  flow  is  guilty  of  generating  as  much  entropy  as  the  entire  entrance 
region?  The  answer  to  this  question  is  simplified  substantially  if,  in  addition 
to  using  Sparrow's  analytical  solution,  we  focus  on  the  t  -*■  0  limit  in  which  the 
viscous  effects  totally  dominate  the  rate  of  entropy  generation  in  the  channel 
(see  the  discussion  following  equation  (2.4)).  Consequently,  to  calculate  S'" 
we  need  only  the  velocity  solution  which  is 

v  =  V  V  *(2i  -  t2)  (2.47) 

X  o  1 


with 

r.  =  y/fi  ,  6/a  =  3(1  -  v^-1) 


(2.48,49) 
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V*  =  V  /V  ,  x* 
1  1  o 


(2.50,51) 


As  shown  on  Fig.  2.5,  V  is  the  uniform  fluid  velocity  at  x  =  0  ,  while  V 

o  1 

anti  6  are  the  centerline  velocity  and  boundary  layer  thickness,  respectively. 
For  the  function  V^*(x*),  Sparrow  reports  a  differential  equation  which, 
integrated,  yields 


x*  =  _L  ( 

x  10  \ 


9V  *  -  2 - —  -  16  In  V, 

1  1 


(2.52) 


Tin?  entrance  region  extends  to  x*  =  0.1038,  corresponding  to  V^*  =  3/2  as 
for  ['lane  Poiseuille  flow. 

Using  equation  (2.1)  and  leaving  out  the  conductive  contribution  we  obtain 


4u  V2 

s'"  =  v^a-u2 


(2.53) 


and,  integrating  across  the  duct. 


2  I  S'"  dy  =  —  — —  V  *2 

J  b  Y  3  T<5  1 


(2.54) 


From  (2.54)  we  find  that  in  the  fully-developed  section  of  the  duct  (FD)  the 
rate  of  entropy  generation  per  unit  area  of  duct  wall  is 


S"  =  6  - 

FD  Ta 


(2.55) 


Fig.  2.')  shows  the  variation  of  S"/S^d  with  axial  position  along  the  duct. 

As  expected,  the  irreversibility  effects  are  most  intense  near  x  =  0  .  However, 


as  x*  increases,  S"  rapidly  approaches  the  fully-developed  entropy  generation 
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flat  rectangular  duct  in  laminar  flow. 


level  S"  a  good  distance  before  the  end  of  the  entrance  region,  x*  =  0.1038. 
FD 

The  total  irreversibility  associated  with  entrance  effects  is  estimated 
integrating  expression  (2.54)  from  x*  =  0  to  x*  =  0.1038  .  The  result  of 
this  operation  is 


11  pa  V 


(2.56) 


where  p  is  the  fluid  density.  Similarly,  the  total  rate  of  entropy  generation 

over  a  fully-developed  stretch  of  length  L  is 

FD 


pv  • 

SFD  =  6  Ta .  LFD  • 


(2.57) 


Comparing  results  (2.56)  and  (2.57),  we  conclude  that  the  fully-developed  equivalent 
of  the  entire  entrance  region  is  a  fully-developed  section  of  length  LfD  given 


Re  =  0.00764  Re  , 
1440  D  D 


(2.58) 


where  D  is  the  hydraulic  diameter  4a  ,  while  ReD  =  VqD/v  •  It  is  time  now 

to  compare  the  irreversibility-equivalent  length  L  with  the  physical  extent 

FD 

of  tile  entrance  region,  L  .  Setting  x*  =  0.1038  and  x  =  L  in  expression 

E  E 

(2.51)  we  find  that  familiar  result 


—  -=  0.00649  Re  . 
D  D 


(2.59) 


Therefore,  from  (2.58,59),  the  irreversibility  contributed  by  the  entrance 
region  is  roughly  equal  to  (actually,  only  about  18%  higher  than)  the  irreversibility 
estimated  assuming  fully-developed  flow  over  a  duct  length  equal  to  the  entrance 
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length.  This  conclusion  is  important  in  practical  irreversibility  calculations 
as  it  is  considerably  easier  to  treat  the  entire  duct  as  in  the  fully-developed 
regime. 
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Nomenclature 


a 


B  ,  B 
o 


c 

V 

C. 

D 


C 


f  ,X 


D 


Ec 

f 


F 

D 

h 

x 

k 


FD 


V  '  NS"'NS' 


NU 


P 

Pr 

q"  .  q'  .  q 

r 


r 

o 

1< 

r<o 


half-thickness  of  flat  duct 
duty  parameters 

specific  heat  at  constant  pressure 

drag  coefficient 

local  skin  friction  coefficient 

hydraulic  diameter 

Eckert  number 

function,  equation  (2.26);  friction  factor,  equation  (2.14) 
drag  force 

local  heat  transfer  coefficient 
thermal  conductivity 
entrance  region  length 

length  of  irreversibility-equivalent  fully  developed  section 

mass  flow  rate 

entropy  generation  number 

Nusselt  number 

pressure 

Prandtl  number 

heat  transfer  interaction,  [W/m2],  [W/m] ,  [W] 
radial  position 
tube  radius 

dimensionless  radial  position 


S,B 


,  S",  S' 


Reynolds  number 
specific  entropy 

rate  of  entropy  generation,  [W/m3K] ,  [W/m2K] ,  [W/mK] 
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S" 

FD 


L 


T 


T 

o 


v  ,  v 
x  y 


V 


y 

a 

6 

<5 

r 

n 

o 

0( 

Vi 

V 


p 

T 


T 

o 


<J> 


rate  of  entropy  generation  in  the  fully-developed  region 
time 

absolute  temperature 

reference  temperature 

specific  internal  energy 

velocity  components 

entrance  velocity 

centerline  velocity 

dimensionless  centerline  velocity 

horizontal  coordinate 

dimensionless  coordinate 

vertical  coordinate 

thermal  diffusivity 

velocity  boundary  layer  thickness 

thermal  boundary  layer  thickness 

dimensionless  coordinate  across  flat  duct 

similarity  variable  in  boundary  layer  flow  over  flat  plate 
temperature  difference,  T  -  T 

o 

extreme  temperature  difference, 
viscosity 

kinematic  viscosity 
fluid  density 

ratio  of  characteristic  temperature  difference  divided  by 
the  absolute  temperature 

wall  shear  stress 

viscous  dissipation  function 
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3.  Till-',  IMPACT  OP  HEAT  TRANSFER  AUGMENTAT  I  ON  ON  ENTROPY  GENERATION 


3 . 1  Entropy  Generation  Analysis 

Consider  a  heat-exchanger  passage  of  length  dx  ,  heat  transfer  per  unit 
length  q’  ,  and  mass  flowrate  m  .  The  passage  geometry  is  described  in  terms 
of  the  hydraulic  diameter  D  and  flow  cross  sectional  area  A  .  The  rate  of 
entropy  generation  per  unit  length  is  [3.l] 


S' 


q*  AT  in  /  dP  \ 
T2  +  pT  \  dx  / 


(3.1) 


where  the  first  term  represents  the  irreversibility  rate  due  to  heat  transfer 
across  the  wall-fluid  temperature  difference,  S '  ,  while  the  second  term  is 

the  irreversibility  rate  attributed  to  fluid  friction,  .  Implicit  in 

writing  Eq.  (3.1)  is  the  assumption  that  the  wall-fluid  AT  is  considerably 
smaller  than  the  local  (wall,  fluid)  absolute  temperature.  This  is  a  good 
assumption,  especially  in  heat  exchangers  operating  above  room  temperature. 

An  important  dimensionless  parameter  in  the  second  law  analysis  ol  con¬ 
vective  heat  transfer  is  the  irreversibil ity  distr ibut ion  ratio 


$  = 


•  /g 

ap  At  ' 


(3.2) 


The  parameter  41  describes  the  relative  importance  of  fluid  friction  in  the 
total  irreversibility  of  a  flow  passage.  Like  and  S^T  ,  the  irreversi¬ 

bility  distribution  ratio  $  is  a  local  parameter. 

The  impact  of  an  augmentation  technique  on  the  irreversibility  of  a  qiven 
heat  exchanger  passage  can  be  evaluated  by  calculating  the  entropy  generation 
rate  in  the  "augmented"  passage,  ,  and  comparing  it  with  the  entropy 
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gonciution  rate  in  the  original  (unaugmented)  passage,  S'  .  It  is  convenient 

o 

to  define  the  augmentation  entropy  generation  number  as  the  ratio 

N  =  S'/S*  .  0.3) 

b ,  a  a  o 

Augmentation  techniques  yieldinq  values  of  N.  less  than  unity  are  thermo- 

S ,  a 

dynamically  advantageous  since,  in  addition  to  enhancing  heat  transfer,  they 
reduce  the  rate  of  entropy  qenerat ion  in  the  apparatus.*  It  is  easy  to  show 
that  when  q'  and  m  are  the  same  before  and  after  augmentation,  the  entropy 
generation  number  can  be  rewritten  as 


N. 

S  ,  ci 


(N  +  <j>oNp)/(l  +  <t>o)  , 


(3.4) 


where 


N  =  (St  /St  ) (D  /D  ) 
T  o  a  a  o 


(3.5) 


N  =■-  (f  /f  )  (D  /D  )  (A  /A  )  2  ,  0.6) 

p  a  o  o  a  o  a 


Re  =  (D  A  /D  /A  )  Re  . 
a  a  o  o  a  o 


(3.7) 


In  these  expressions,  St  ,  f  ,  Re  ,  D  and  A  are 
fact or,  and  Reynolds  number,  hydraulic  diameter  and 


the  Stanton  number,  friction 
flow  cross  sectional  area, 


*  lor  example,  in  the  regenerative  heat  exchanger  for  a  Brayton  cycle  heat 
engine  a  low  Ng>a  means  less  dissipation  of  available  work  in  the  power  cycle. 
Among  the  physical  parameters  affected  directly  are  the  stream-to-stream  tempera¬ 
ture  difference  and  the  two  pressure  drops  on  either  side  of  the  heat  transfer 
surface.  As  the  heat  exchanger  entropy  generation  decreases,  the  pressure  drop 
icross  the  turbine  increases,  resembling  more  closely  the  pressure  rise  provided 
by  the  compressor.  Consequently,  the  turbine  power  output  increases.  At  the 
same  time,  the  inlet  temperature  to  the  combustion  chamber  (heater)  increases, 
hence,  the  specific  fuel  consumption  decreases.  We  see  that  the  reduction  in 
entropy  generation  in  the  counterflow  heat  exchanger  (reduction  indicated  by 

Nr.  -1)  can  be  measured  directly  as  an  increase  in  energy  conversion  efficiency. 
,  a 
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* 

respectively.  Equation  (3.4)  shows  that,  when  the  reference  passage  is 

dominated  by  heat  transfer  irreversibility  0)  ,  the  augmentation 

number  N  reduces  to  N  which  is  proportional  to  the  ratio  of  heat  trans- 
S  §  &  T 

fer  coefficients  (reference/augmented).  In  general,  however,  the  augmentation 

entropy  generation  number  N  is  a  function  of  the  ratio  of  heat  transfer 

S  f  ct 

coefficients  and  the  friction  factor  ratio. 

In  what  follows,  we  use  Eq.s  (3.4)-(3.6)  to  evaluate  the  relative  change 
in  irreversibility  introduced  by  each  of  the  augmentation  techniques  of 
Fiq.  3.1.  We  obtain  this  result  by  using  published  heat  transfer  (St)  and 
fluid  friction  (f)  information  on  each  of  these  flow  geometries. 

5 . 2  Twisted  Tape  Inserts 

The  use  of  twisted  tape  inserts  to  auqment  both  laminar  and  turbulent  con¬ 
vective  heat  transfer  is  one  of  the  most  common  in-tube  augmentation  techniques. 
This  technique  is  also  well-documented.  For  example,  Lopina  and  Bergles  [3.2] 
measured  the  heat  transfer  and  fluid-friction  characteristics  of  a  tightly 
inserted  uninsulated  tape  in  turbulent  flow.  The  tape  was  inserted  in  a  straight 
tube  ami  the  assembly  was  redrawn,  with  tape  in  place.  This  operation  forced 
the  tape  to  slightly  penetrate  the  interior  surface  of  the  tube.  The  geometry 
of  the  augmented  arrangement  is  shown  schematically  in  Fig.  3.1(a).  The 
experimental  results  for  St  and  f  agree  to  within  10%  with  the  published 

cl  c3 

work  of  other  experimenters. 

Figure  3.2  illustrates  the  augmentation  entropy  generation  number  for 
twisted  tape  inserts  in  turbulent  flow,  based  on  the  data  of  Lopina  and  Bergles 
[:..!].  On  this  basis,  we  obtain 

y  =  2.45,  Nt  =  0. 49/(1  +  0.89  Re^0,13)  ,  Np  =  2.33  ;  (3.8) 
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y  =  5.08 


0. 44/(1  +  1.3  Re 


-0.13 


Np  =  3.13 


(3.9) 


also 


D  /D  =  0.61  ,  A  /A  =  1  . 
a  o  a  o 


In  this  example,  the  geometry  of  Fig.  3.1(a)  is  the  augmented  design,  vis-a- 
vis  the  smooth  straight  tube  which  is  used  as  reference. 

It  is  evident  from  Fig.  3.2  that  the  parameter  which  influences  N 

S ,  a 

the  most  is  the  irreversibility  distribution  ratio  4>  .  Thus,  whenever  4* 

o  o 

is  less  than  approx.  0.25,  the  use  of  twisted  tape  inserts  will  lead  to  a 

reduction  in  thermodynamic  irreversibility.  The  greatest  reduction  occurs 

when  4>  =  0  ,  i.e.,  when  the  fluid  friction  contribution  to  S'  is  totally 

o 

negligible  and  N  becomes  proportional  to  the  ratio  of  heat  transfer 
S  /  cl 

coefficients,  [see  Eq.  (3.5)].  It  is  important  to  note  that  the  critical  value 
of  4>o  above  which  the  augmentation  technique  will  actually  increase  the  rate 
of  entropy  generation,  4>q  -  0.25  ,  corresponds  to  a  heat  exchanger  passage 
in  which  the  fluid  friction  contribution  amounts  to  only  20%  of  S'  .  There¬ 
fore,  even  in  cases  clearly  dominated  by  heat  transfer  irreversibility  (4>  <  1) 

o 

it  is  important  to  evaluate  the  total  entropy  generation  rate  and  the  change 

induced  in  it  by  the  implementation  of  a  proposed  augmentation  technique, 

before  deciding  whether  heat  transfer  augmentation  is  thermodynamically  advantageous. 

Figure  3.2  shows  also  that  other  parameters,  such  as  Re^  and  y  ,  are 

not  nearly  as  important  in  affecting  the  value  of  N  .  In  particular,  the 

b  (d 

augmentation  entropy  generation  number  is  almost  independent  of  the  reference 

Reynolds  number,  Re^  .  This  behaviour  follows  from  the  fact  that  the  measured 

Re  dependence  of  St  and  f  if  functionally  the  same  as  that  of  St 

a  a  a  o 

and  f  (see  Refs.  [3.2]  and  [3.3]). 
o 
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Thu  twist  ratio  y  has  a  visible  effect  on  the  entropy  generation  rate; 

for  example,  when  approaches  zero  a  tighter  twist  (y  =  2.45)  is  preferable. 

On  the  other  hand,  for  values  of  p  other  than  in  the  $  ->  0  limit,  i.e., 

o  o 

when  the  fluid  friction  irreversibility  begins  to  play  a  role,  a  looser  twist 
is  more  advantageous  (y  =  5.08).  In  i act ,  using  a  more  extensive  collection 
of  experimental  data,  it  may  be  possible  to  determine  the  required  twist  ratio 
y  which,  for  a  given  ,  is  capable  of  lowering  ^  below  a  predetermined 

value. 

3.3  1  lelical  Tubes 

Another  heat  transfer  augmentation  feature  common  in  advanced  heat  exchanger 
design  is  the  use  of  spiraled  tubes.  In  this  configuration  the  swirl  flow 
induced  in  the  stream  filling  the  tube  is  the  result  of  centrifugal  forces. 

As  shown  iii  Fig.  3.1(b) ,  the  geometry  is  described  by  the  ratio  X  =  R/r  . 

Figure  3.3  was  constructed  by  using  the  heat-transfer  and  fluid-friction 
measurements  reported  by  Seban  and  McLaughlin  [3.4], 

N  =  (X2/Re  )0’05  ,  N  =  (Re  /X2)0’05  ,  D  /D  =  1  ,  A  /A  =  1  . 

T  o  Po  ao  ao 

(3.10) 

in  Fig.  3.3,  we  show  the  augmentation  entropy  generation  number  in  turbulent 
flow  for  only  two  geometries,  X  =  17  and  X  =  104  .  For  brevity,  the 
analysis  leading  from  the  Seban  and  McLaughlin  correlations  [3.4]  to  Eq.  (3.10) 
is  not  shown.  However,  the  reader  can  find  an  outline  of  this  analysis  in 
Ref.  [3.3]. 

The  chief  conclusion  to  be  drawn  from  Fig.  3.3  is  that  shaping  the  tube 
as  a  helix  has  only  a  minor  effect  on  the  rate  of  entropy  generation  in  the 
tube.  In  the  limit  -*  0  ,  the  irreversibility  is  reduced  to  roughly  50%  of 
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the  original  (reference)  value.  In  the  opposite  extreme,  when  fluid  friction 

dominates,  (4>  -*■  >»)  ,  entropy  generation  number  is  of  order  one.  This  means 

o 

that  the  coiling  of  the  tube  has  little  effect  on  the  overall  pressure  drop 

across  the  tube.  Consequently,  for  a  wide  range  of  values  of  $  the  use 

o 

of  helical  tubes  leads  to  a  reduction  in  entropy  generation. 

3. 4  Propeller  Inserts 

A  swirl  flow  promoting  technique  similar  to  twisted  tape  inserts  is  shown 
schematically  in  Fig.  3.1(c).  This  augmentation  technique  consists  of  a  sequence 
of  propellers  placed  rigidly  inside  a  circular  duct. 

In  this  section  we  evaluate  the  augmentation  entropy  generation  number 
corresponding  to  the  air  experiments  reported  by  Koch  [3.5].  The  propeller¬ 
shaped  inserts  investigated  by  Koch  were  fabricated  by  making  six  radial  cuts 
(almost  to  the  center)  in  a  50  mm  diameter  disc.  The  resulting  blades  were 
twisted  to  approx,  a  45°  angle.  A  brass  rod  was  then  inserted  through  the 
center  of  each  disc  and  the  rod  and  discs  assembly  was  soldered,  maintaining 
a  disc-to-disc  spacing  L'  .  The  solid  assembly  was  then  inserted  in  a  straight 
tube  and  secured  in  place  to  prevent  rotation. 

Using  Koch's  data  [3.5]  we  derive  the  following  extreme  values  of  the 
augmentation  entropy  generation  number,  Nt  and  Np  : 

L'  =  49  mm  ,  N  =  0.035  Re  °'15  ,  (L'/D  =  1)N  =  59.8  Re  _0'009  ,  d  /D  =  0.66 

T  o  op  o  a  o 

(3.11) 

L'  =  98  mm,  N  =  0.078  Re  °'114  (L'/D  =  2)N  =  32.4  Re  -0-011  D  /D  =  0.80  ; 

T  o  o  P  o  a  o 

(3.12) 

L'  =  325  mm  ,  N  =  0.153  Re  °'085  ,  (L'/D  =  6.5)N  =  8.6  Re  0 ‘ 094  ,  d  /D  =  0.93 


(3.13) 
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L*  =  980  nun  ,  N  =  0.21  Re  °-085  (L'/D  =  19.6)N  =  3.83  Re  0-11  ,  D  /D  =  0.98  ; 

T  o  O  P  o  a  o 


(3.14) 


where  A  =  A  . 

a  o 

Figure  3.4  shows  the  augmentation  entropy  generation  number  for  three  values 

of  irreversibility  distribution  ratio,  <j>  =  0,  1,  °°.  As  in  the  preceding  examples, 

the  dependence  of  N  on  Re  is  weak.  However,  the  maximum  reduction  in 

S  ,a  o 

entropy  generation  is  more  substantial  than  in  the  first  two  examples.  For  a 

relatively  dense  sequence  of  propellers,  L'  =  49  mm  ,  in  the  limit  d>  ■+  0  the 

o 

entropy  generation  rate  is  reduced  by  as  much  as  a  factor  of  ten.  Figure  3.4  also 
shows  that,  depending  on  L'  ,  the  reference  irreversibility  distribution  ratio 
<j>  must  be  less  than  0.05  if  propeller  inserts  are  to  be  effective  in  conserving 
available  work. 

3. 5  Internally  Finned  Tubes 

Figure  3.1(d)  shows  schematically  a  round  tube  covered  internally  with  straight 
or  spiraled  fins.  This  augmented  tube  geometry  has  been  studied  extensively.  For 
example,  Bergles  et  al.  [3.6]  and  Carnavos  [3.7]  have  measured  the  heat  transfer 
and  fluid  friction  characteristics  of  approximately  thirty  straight  and  spiraled 
internally  finned  geometries.  Empirical  formulae  generated  by  Carnavos  [3.7] 
correlate  within  10%  the  bulk  of  the  existing  experimental  information. 

In  this  section,  we  use  the  Carnavos  correlations  to  examine  the  relationship 
between  internal  finning  and  entropy  generation.  A  first  objective  is  to  determine 
the  difference  between  straight  and  spiraled  fins  of  otherwise  similar  geometry. 

The  effect  of  fin  density  and  twist  ratio  is  also  illustrated. 

Figure  3.5  shows  the  change  in  entropy  generation  rate  due  to  straight  internal 
fins  in  turbulent  flow,  for  the  two  extremes  <j>  =  0,  °°  .  Three  different  geometries 


wamvmmmm 
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have  been  considered,  n  =  6,  10,  12,  where  n  is  the  number  of  fins.  The  re¬ 
sults  of  the  entropy  generation  analysis  which  led  to  Fig.  3.5  can  be  summarized 
as  follows  [3.3] : 


for 

3 

O' 

N  =  0.097 

T 

„  0.2 
R6o  ' 

N  =  7.35  Re  °'1  , 
P  0 

D  /D  =  0.74  ,  A  /A  =  0.94  ; 
a  0  a  0 

(3.15) 

for 

n  =  10  , 

N  =  0.089 
T 

Re  °-2 

O 

,  Np  =  9.45  Reo~°* 1  , 

D  /D  =  0.61  ,  A  /A  =0.91 
a  0  a  0 

(3.16) 

for 

n  =  12, 

N  =  0.084 

T 

0.2 

Reo  ' 

N  =  10.6  Re  _0'1  , 

P  0 

D  /D  =  0.57  ,  A  /A  =  0.89  . 
a  0  a  0 

(3.17) 

It  is  evident  that,  almost  independently  of  n  ,  the  reduction  in  entropy  genera¬ 
tion  is  small  even  in  the  limit  4>q  -*■  0  ,  particularly  as  the  Reynolds  number 
increases .  The  change  in  entropy  generation  rate  is  more  pronounced  as  the  number 
of  fins  increases. 

Greater  reductions  in  entropy  generation  are  achieved  when  the  internal  fins 
are  twisted,  thereby  inducing  a  swirl  motion  in  the  core  of  the  tube  cross  section. 

In  Fig.  3.6,  we  show  the  augmentation  entropy  generation  number  corresponding  to  a 
finned  tube  with  n  =  10  and  the  values  of  twist  ratio,  y  =  5.62  and  y  =  9.27. 

In  general,  the  augmentation  entropy  generation  number  is  given  by  Eq.  (3.4) 
combined  with 

N  =  0.052  Re  0.17  0.27  _  Re  ,  D  /D  =  0.61  ,  A  /A  =  0.91  . 

T  o  P  o  a  o  a  o 

(3.18) 

Provided  the  irreversibility  distribution  ratio  i)>  is  small  enough,  the  reduction 
in  entropy  generation  rate  increases  as  y  decreases,  i.e.,  as  the  fin  twist 


Fig.  3.6. 


’ig.  3.7. 


Augmentation  entropy  generation  number  for  internal  spiraled  fins. 


10 


1.0 


0.1 

i04  Id*  R)So  IC f 

Comparison  of  four  swirl  promoters  based  on  irreversibility  minimiza 
tion  potential  in  a  <{>^=0.1  application;  IF=internal  straight  fins 

(n~10) ;  lSF=internal  spiraled  fins  (n=10,y=9.27) ;  TT=twisted  tape 
(y-5.08);  H=helical  tube  (A=30) . 


-41- 


becomes  tighter.  This  conclusion  is  qualitatively  similar  to  the  one  encountered 

earlier  in  the  analysis  of  twisted  tape  inserts. 

Comparing  the  curves  of  Fig.  3.6  with  the  n  =  10  curves  of  Fig.  3.5,  we 

conclude  that  the  twist  of  internal  fins  can  have  a  significant  effect  on  the 

available  work  conservation  potential  of  the  finning  technique.  For  example, 

in  the  limit  <|>  -*•  0  ,  the  spiraled  fin  irreversibility  is  about  50%  less  than 

o 

in  the  case  of  straight  fins. 

An  imjxirtant  observation  is  that  the  spiraled  fins  entropy-generation  number 

N  (<f>  ,  Re  ,  y)  ,  given  by  Eqs.  (3.4)  and  (3.18),  is  not  a  monotonic  function 
S,a  o  o  ' 

of  twist  ratio  y  .  In  fact,  when  the  heat  exchanger  passage  is  specified  (i.e., 

a  and  Ro  are  fixed) ,  one  can  determine  the  optimum  twist  ratio  y 
o  o  opt 

yielding  the  largest  reduction  in  entropy  generation  rate  (minimum  N  )  .  Setting 

S ,  a 

IN  /3y  =  0  and  3N  /3y  >  0  ,  we  obtain 

•’  /  S  f  3 


„  .  ,  2.13  n  -0.64 

y  =  9.4  x  104  ■p  Re 

opt  o  o 


(3.19) 


For  example,  in  an  application  in  which  $  =0.2  and  Re  =  104  ,  the  optimum 

o  o 

twist  ratio  is  y  =  8.4  ,  yielding  the  minimum  N™10  =  0.868.  If,  in  the 
opt  1  S , a 

same  application,  heat  transfer  is  augmented  using  fins  with  twist  ratios 


y  /  y 


opt 


say  y  =  5  or  y  =  15  ,  the  entropy  generation  reduction  is  inferior 


and  N  =  0.875  and  N  =  0.876  ,  respectively. 

S  /  3  S  $  cl 


3 . 6  Internally  Roughened  Tubes 

Finally,  we  illustrate  the  evaluation  procedure  in  the  case  of  "in-tube 
roughness)  which  is  one  of  the  most  common  and  amply  documented  heat  transfer 
augmentation  techniques.  In  this  case  the  augmentation  entropy  generation  number 


(3.4)  assumes  the  simpler  form 
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N 


S  f  ci 


St  0 

— 2.  +  - 2 — 

St  1  +  0 
a  o 


(3.20) 


since,  in  almost  all  rough-tube  designs,  and  A^  -  A^  .  Consequently, 

the  Reynolds  number  in  the  augmented  (rough)  tube  is  equal  to  the  Reynolds 
number  in  the  reference  (smooth)  tube.  The  entropy  generation  number  for  in¬ 
tube  roughness  emerges  as  a  weighted  function  of  the  Stanton  number  ratio  and 
the  friction  factor  ratio.  The  two  ratios  are  weighted  according  to  the  irrever¬ 
sibility  distribution  parameter  0  ,  in  other  words,  with  respect  to  the 

o 

irreversible  mode  of  operation  of  the  heat  exchanger  duct  (fluid  friction  vs. 
heat  transfer) . 

Figure  3.9  shows  the  entropy  generation  number  as  a  function  of  Re  and 

o 

0  .  The  figure  was  constructed  for  Pr  =  0.72  using  the  data  of  Refs.  [3. 8, 3. 9] 

o 

for  sand  grain  (f  ,  St  )  and  the  correlation  of  Refs.  [3.10,3.11]  for  smooth 

tubes  (f  ,  St  ) .  It  is  evident  that  the  ability  of  sand-grain  roughened  walls 
o  o 

to  reduce  irreversibility  depends  on  0q  and  Re^  .  The  effect  of  roughness 

height  e/D  on  entropy  generation  increases  steadily  as  the  Reynolds  number 

increases.  As  one  might  expect,  the  change  in  N  is  more  pronounced  in  the 

S ,  a 

case  of  larger  e/D  . 

From  an  engineering  standpoint,  it  is  important  to  know  under  what  conditions 

the  wall- roughening  leads  to  a  reduction  in  entropy  generation.  In  Fig.  3.10 

we  plotted  the  marginal  curves  0  (Re  )  resulting  from  setting  N  =  1  in 

o  o  S ,  a 

equation  (3.20).  Below  each  curve,  roughening  the  wall  will  decrease  the  irrever¬ 
sibility  rate.  Above  each  curve,  wall-roughening  will  increase  the  rate  of 
entropy  generation. 

Figures  3.11-3.14  show  the  augmentation  entropy  generation  number  N  and 

S  1 3 

marginal  curves  0  (Re  )  for  repeated  ribs  of  the  type  sketched  in  the  insert 
o  o 

of  Fig.  3.11.  This  geometry  is  characterized  by  two  aspect  ratios,  the  rib 


-*r 


r~  “  r  r  i  i  t  i  t  17  -  -  t  r — r  t  r  •  rq  1  1  »  '  '  '  '  1 


\<f  10*  10P  \<f  I  o' 


Augmentation  entropy  generation  number  for  rind-grain  roughness 


Augmentation  entropy  generation  number  for  repeated  rib  roughnes 
effect  of  rib  height  e/D  . 


Critical  irreversibility  distribution  ratio  <j>  for  repeated 
rib  roughness;  effect  of  rib  height  e/D  .  1 


Augmentation  entropy  generation  number  for  repeated  rib  roughness 
effect  of  rib  spacing  P/e  . 


Critical  irreversibility  distribution  ratio  $  for  repeated 
nb  roughness;  effect  of  rib  spacing  P/e  .  ° 
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height  e/D  and  the  rib  spacing  P/e  .  The  graphs  were  constructed  for 

Pr  ■-  0.72  using  the  friction  factor  and  Stanton  number  correlations  of  Webb, 

e t  a  1  ■  [3.12].  There  are  obvious  similarities  between  these  plots  and  Figs.  3.9 

and  3.10  for  sand-grain  roughness.  The  more  significant  difference  is  that, 

with  Ke^  and  e/D  fixed.  Figs.  3.12  and  3.14  yield  marginal  irreversibility 

distribution  ratios  <j>  about  four  times  less  than  Fig.  3.10.  In  other  words, 

o 

repeated  ribs  of  a  given  height  (e/D)  correspond  to  sand-grain  roughness  with 

grain  height  significantly  greater  than  (e/D)  . 

The  effect  of  rib  spacing  (P/e)  is  summarized  in  Figs.  3.13,  3.14.  We 

f ind  that  the  rib  spacing  has  a  relatively  weaker  effect  on  entropy  generation 

than  the  rib  height  (Figs.  3.11,  3.12).  Furthermore,  the  change  induced  in  the 

entropy  generation  number  decreases  steadily  as  the  rib  spacing  increases.  This 

trend  is  in  agreement  with  the  effect  of  decreasing  e/D  for  sand  grain 

(Fig.  3.9);  increasing  the  rib  spacing  while  keeping  the  rib  height  constant  is 

equivalent  to  gradually  smoothing  out  the  tube  surface. 

An  important  observation  concerns  the  minima  in  the  Ng  ^  curves,  Figs.  3.11, 

1.13,  particularly  in  the  range  in  which  the  use  of  repeated  ribs  leads  to 

a  reduction  in  irreversibility.  For  a  given  heat  exchanger  duct  (4>  ,  Re  )  it  is 

o  o 

possible  to  determine  a  priori  the  optimum  rib  geometry  which  yields  the  maximum 
reduction  in  the  rate  of  entropy  generation.  For  example,  in  Fig.  3.15  we 


report  the  optimum  rib  height  (e/D)  for  cases  where  the  relative  spacing  is 


fixed,  !•/<.’  -  10  .  Note  that  the  irreversibility  distribution  ratio  4>  has  a 

o 

marked  effect  on  the  optimum  rib  height  for  minimum  N  ,  despite  the  fact 

b  /  3 

t hat  numerically  i  is  very  small. 

o 


(en.  ration  m  the  case  of  a  limited  number  of  augmentation  techniques. 


Fig.  3.15.  Optimum  rib  height  for  minimum  exergy  destruction. 

For  a  considerably  more  complete  study  the  reader  should  consult  the  thesis 
by  Ouellette  [3.3] .  In  addition  to  the  techniques  highlighted  in  this 
chapter,  Ouellette's  thesis  analyzes  the  entropy  generation  due  to 

I  -  extended  surfaces 

-  internal  fins 

-  mesh  inserts 

II  -  rouqher.ed  surfaces 

-  thread- type  roughness 

-  screen  roughness 

-  wire-wound  roughness 

III  -  displaced  enhancement  devices 

-  suspended  discs 

-  suspended  rings 
IV  -  compound  techniques 

-  twisted  tape  with  fin  effect 

-  internally  spiraled  fins 


The  analytical  results  of  Ouellette's  study  are  summarized  here  in  Tables  3. 1-3. 5. 
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4.  EXTENDED  SURFACES  (FTNS)  FOR  MINIMUM  ENTROPY  GENERATION 

4 . 1  Entropy  generation  due  to  heat  transfer  from  a  single  fin 

The  entropy  generation  associated  with  a  fin  in  cross-flow  can  be  evaluated 
based  on  the  general  model  presented  in  Fig.  4.1.  Consider,  then,  a  solid  body 
of  arbitrary  shape  and  surface  area  A  ,  suspended  in  a  uniform  stream  of  free- 
streom  velocity  Ua  and  temperature  T^  .  Imagine  also  a  stream  tube  (control 
surface)  surrounding  the  solid  body,  the  radius  of  this  tube  being  considerably 
larger  than  the  characteristic  linear  dimension  of  the  body.  Since  the  fluid 
outside  and  immediately  inside  the  tube  surface  belongs  to  the  free  stream  , 

T  ,  the  stream  tube  surface  can  be  regarded  as  adiabatic  and  shear-free. 

The  equations  describing  the  conservation  of  mass,  conservation  of  energy 
and  generation  of  entropy  in  the  stream  tube  are 


r 

V 
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rate  can  finally  be  written  as 


gen 


l"  (  — - —)  do  +  —  F  u  . 

i  I  T  T  ,  /  T  D  00 

'  OO  W  oo 


(4.5) 


A 


In  order  to  calculate  the  total  rate  of  entropy  generation  due  to  a  single 
fin,  it  is  necessary  to  recognize  that  expression  (4.5)  sums  up  only  the  entropy 
generated  in  the  space  surrounding  the  fin  surface.  It  is  easy  to  see  that 
since  the  is  not  an  isothermal  body,  entropy  will  also  be  generated  internally. 
In  accordance  with  the  second  law  of  thermodynamics  applied  to  the  fin  shown 
in  the  bottom  half  of  Fig.  4.1,  the  entropy  generated  inside  a  single  fin  is 


(S 


gen  internal 


q"  ,  B 
do  -  — 


(4.6) 


wiiere  q  and  T  represent  the  base  heat  transfer  and  absolute  temperature, 
fl  B 

The  total  rate  of  entropy  generation  associated  with  a  single  fin  is  obtained 
by  adding  equations  (4.5)  and  (4.6)  side-by-side. 


q  0  , 

(S  )  =  -S-7—  +  -  F  U 

gCn  fin  T°°  T»  °  “ 


(4.7) 


[n  this  expression  0^  is  the  base-stream  temperature  difference  (T  -  T^)  , 

which  is  assumed  considerably  smaller  than  the  local  absolute  temperature, 

T  or  T  . 

B 

The  fin  entropy  generation  rate,  Equation  (4.7),  is  a  remarkably  simple 

result  which  demonstrates  that  fluid  friction  and  inadequate  thermal  contact 

contribute  hand-in-hand  to  the  degrading  of  fin  thermodynamic  performance. 

The  heat  transfer  contribution  to  (S  )  is  proportional  to  the  product 

gen 

fin 

of  base  heat  flux  times  base-fluid  temperature  difference.  Therefore,  the 
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heat  transfer  entropy  generation  term  can  be  minimized  by  reducing  0  ,  r.e., 

B 

by  enhancinq  the  base-stream  thermal  contact.  In  general,  this  operation 
requires  an  increase  in  fin  size  (surface  area)  which  leads  automatically  to 
an  increase  in  drag  and  fluid  friction  entropy  generation.  Clearly,  the  fin 
size  plays  an  important  trade-off  role  with  respect  to  improving  the  thermo¬ 
dynamic  performance  of  the  individual  fin.  In  what  follows  we  illustrate  this 
trade-off  by  showing  concrete  results  [  4 . 1 ]  for  entropy  generation  minimization 
in  some  of  the  most  common  fins  encountered  in  practice. 

4.2  pin  Fins 


Consider  first  the  pin  fin  geometry  shown  in  Fig.  4.2(a).  This  geometry 
is  one  of  the  simplest,  because  it  depends  on  only  two  dimensions:  the  length 
L  and  the  diameter  of  the  circular  cross-section  D  .  According  to  the  uni¬ 
directional  heat  conduction  model  described  in  the  Introduction,  the  relationship 
between  base  heat  flux  and  base-stream  temperature  difference  is  [4.2] 


TT  o 

—  kD^m  tanh (mL) 

4 


(4.8) 


Substituting  this  expression  into  Equation  (4.7),  the  total  entropy  generation 
rate  can  be  written  as 


p  v2  ReL  Reo  CD 

7T  9  S  S  V 

_  rn  Z  /  \|#\  v  r . *  D  _  f  anV. 

2NU*5  f  A 

Vs  R6l' 

2  T 

00 

2  Too  (Ak)  Nu  ReD  u  tanh 

^  oo 

2NU  ^ 

■j  ReD_ 

(4.9) 


where  the  drag  coefficient  is 


V(?pU~  DL) 


Since  the  external  flow 


is  assumed  known,  we  are  using  the  Reynolds  number  as  dimensionless  notation 


ir  common  fin  geometries:  (a)  pin  fin;  (b)  rectangular  plate  fin;  (c)  rectangular  plate  fin 
trapezoidal  profile;  (d)  triangular  plate  fin  of  rectangular  profile. 
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for  the  two  dimensions  of  the  pin  fin,  Re  =  U  D/v  ,  ReT  -  U  L/v  . 

D  "  L  ° 

In  this  study,  we  adopted  the  design  viewpoint  that  the  "30b"  of  the 

individual  fin  is  to  transfer  heat  at  a  known  rate  (q  )  from  the  wall  to  the 

B 

known  stream,  in  the  least  irreversible  manner  possible  (with  minimum  genera¬ 
tion  of  entropy).  Therefore,  the  correct  entropy  generation  number  [4.3]  for 
this  design  problem  is  constructed  as 


N,.  = 


gen 


fin 


q  ‘  U 

k  v  T  2 


(k/A)S 


r  .  .  ,  ,*5  Re  ] 

ii  5 

5  /  A  \  L 

—  Nu  Re  tanh 

2NU  - 

2  D 

\  k  /  Re  , 

L  D J 

Re 


L 


Re 

D 


(4.10) 


where  B  is  a  fixed  dimensionless  parameter  which  accounts  for  the  importance 
of  fluid  friction  irreversibility  relative  to  heat  transfer  irreversibility , 


B  =  p  v  *  k  T  /q  2  . 

O)  ^3 


(4.11) 


It  tne  pin  fin  is  slender,  the  Nusselt  number  and  the  drag  coefficient  can  be 
evaluated  from  results  developed  for  a  single  cylinder  in  cross-flow  [4.4] 


Nu  =  0.683  Re  0-466  pr1,/3  ,  C  =  5.484  Re  0-246  40  <  Re  <  103 

D  D  D  '  D 


(4.12) 


The  entropy  generation  number  N<_  emerges  as  a  function  of  five  dimension- 

loss  groups,  two  pertaining  to  fin  geometry  (Re  ,  Re  )  ,  and  three  accounting 

Li  D 

to r  the  working  fluid  and  for  the  fin-stream  convective  arrangement  (Pr,  k/\,  B) . 

Minimization  of  with  respect  to  Re^  is  achieved  in  a  straightforward 

manner  by  solving  <)N  /3Re  -  0  .  The  optimum  pin  length  calculated  in  this 

o  L 


manner  is 
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Ke 


opt 


Re 
_ D 

i. 

JNu 


s 

r  S  n 

/  k  \  -1 

/  »  \ 

(  -  -  I  s inh 

(  - 1 “ - 1“  ) 

\  A  } 

n  Re  '  t  B  / 

L'  D  D  J 

(4.13) 


The  engineering  sign  1  f icance  of  result  (4.13)  is  that  the  optimum  pin 

length  cun  be  calculated  immediately,  provided  Re^  (hence,  Nu  and  C  )  is 

specified.  In  most  cases,  however,  the  pin  diameter  is  a  design  variable  which 

can  also  be  determined  on  the  basis  of  entropy  generation  minimization.  The 

optimum  diameter  for  minimum  entropy  generation  Re  can  be  determined  by 

opt 

substituting  equations  (4.12)  and  (4.13)  into  equation  (4.10),  and  numerically 


minimizing  N 

S 

An  alternative  approach  to  sizing  a  pin  fin  for  minimum  irreversibility 

consists  of  determining  the  optimum  diameter  Re^  subject  to  fixed  slender- 

opt 

ness  ratio  y  -  Re  /Re  .  This  constraint  stems  from  practical  limitations 

Li  U 

encountered  in  the  process  of  manufacturing  a  surface  covered  with  a  large 

number  of  fins.  The  entropy  generation  number  (4.10)  can  be  expressed  as  a 

function  of  Re^  and  y  .  Representative  results  of  the  numerical  work  of 

minimizing  are  shown  in  Fig.  4.3.  The  entropy  generation  number  has 

u  clear  minimum  with  respect  to  the  pin  diameter  Re^  ,  when  y  ,  B  ,  and  M 
* 

are  fixed.  The  optimum  pin  diameter  Re  increases  if  the  slenderness 

D ,  opt 

ratio  decreases.  This  general  trend  is  summarized  in  the  right  half  of  Fig.  4.3, 

—8  —5 

for  t  lie  range  5  <_  y  15  ,  10  B  10  .  When  the  slenderness  ratio  is 

fixed,  the  optimum  pin  diameter  decreases  as  the  fluid  friction  becomes  more 
important  in  the  entropy  generation  total  (i.e.,  as  B  increases). 

Another  practical  design  constraint  worthy  of  consideration  is  the  fin 
volume.  This  constraint  is  forced  upon  the  designer  by  the  high  cost  of  con¬ 
struction  materials  such  as  copper.  A  dimensionless  group  which  is  proportional 

to  the  volume  of  the  pin  fin  is  V  =  Re  2  Re  .  The  entropy  generation  number 

D  b 


M  -  (k/i)’,/Pr1/6 


A 
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(N^)  can  bo  expressed  in  terms  of  V  and  another  geometric  variable,  for  ex¬ 
ample,  tlie  slenderness  ratio.  on  the  left  side  of  tig.  4.4  we  show  the  genera 
behavior  of  the  entropy  generation  number  ,  as  y  and  B  vary  subject 

to  the  constant  volume  constraint.  The  optimum  slenderness  ratio  y  for 

opt 

minimum  N.  is  summarized  in  the  right  half  of  Fig.  4.4:  above  a  certain 


value  of  friction  parameter  B 
!t  increases. 


the  optimum  ratio  y 

opt 


drops  off  rapidly  as 


1 


4 . 3  Plate  Fins 

Ln  this  section  we  focus  on  fin  geometries  which  can  be  modeled  as  thin 
conducting  plates  parallel  to  the  flow  direction.  Fins  of  this  type  vary  widely 
with  respect  to  the  shape  and  cross-section  of  the  surface  swept  by  the  flow: 
m  many  cases  the  tin  is  wide  at  the  base  and  narrow  at  the  tip,  in  approximate 
agreement  with  Schmidt's  principle  of  material  reduction  in  a  fin  with  fixed 
heat  transfer  [4.5],  below,  we  consider  the  minimum  entropy  generation  design 
of  three  frequently  used  plate  fin  geometries,  illustrated  in  Figs.  4.2(b), 

(c) ,  and  (d) . 

Rectangular  Plate  tin.  The  simplest  plate  fin  geometry  is  shown  schematically 
in  Fig.  2(b).  The  minimum  irreversibility  design  of  this  fin  requires  the 
selection  of  three  geometric  parameters,  the  length  h  ,  the  breadth  (length 
swept  by  fluid)  b  ,  and  the  plate  thickness  <5  ,  where  L  »  b  >>  6.  As  in 
tin  preceding  section,  we  are  placing  the  study  in  the  limit  where  the  plate 
tin  is  lender  (b  •-  I.);  in  this  limit,  we  rely  on  laminar  heat  transfer  and  skin 

friction  results  for  two-dimensional  flat  plates  [4.4]: 
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The  entropy  generation  number  for  a  rectangular  plate  fin  in  laminar  flow  be¬ 
comes 
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1.15  Re.*5  Re,  Pr tanh  (mL) 
o  b 


+  1.328  B  Re  Re, 


(4.15) 


where 


B  =  p  v  - 


k  TyqB* 


mb  =1.15 


(*) 


1/6  -1/4 

Pr  Re  Re,  Re 

L  b 


—  *2 


(4.16) 


From  the  outset  we  note  that  the  plate  thickness  6  appears  only  in  the 
heat  transfer  term  of  ,  consequently,  Re^  does  not  play  a  trade-off  role 
in  the  minimization  of  Ng  .  Since  in  most  practical  applications  6  is 
determined  by  considerations  such  as  price,  availability  and  machinability  of 
sheet  metal,  it  makes  engineering  sense  to  regard  Re^  as  fixed.  The  minimi¬ 
zation  of  N  with  respect  to  Re  and  Re  is  achieved  by  solving  the 
S  b  L 

simultaneous  set  of  equations  3N  /3Re  =  0  and  3N  /3Re  =  0  :  the  optimum 

S  b  S  L 

geometry  for  minimum  entropy  generation  is  explicitly  given  by 


-2/3  -2/3 

Ro  =  0.984  B  Re,  ,  Re 

b,opt  6 


„  -  0.685  (*/«'  R.//3  Pr'VV1/6 

L,opt  0 

(4.17) 


This  also  means  that  the  optimum  slenderness  ratio  y  =  L/b  is 


y  =  Re  /Re,  =  0.696  MB  Re. 
opt  L  b  6 

opt  opt 


(4.18) 


t n  conclusion,  the  optimum  plate  fin  dimensions  b  and  L  can  be 
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calculated  directly  using  equations  (4.17)  as  soon  as  the  base  heat  flux  (q  ), 

6 

the  flow  field  and  the  sheet  metal  thickness  ( <5 )  are  known.  It  is  necessary 
to  keep  in  mind  that  the  present  formulas  are  valid  only  in  the  laminar  regime, 


lu  ,  <  5  x  10  ,  and  in  the  slender  shape  limit  y  ^  1  .  Similar  re- 

n  opt 

opt 

suits  can  be  developed  numerically  for  the  turbulent  regime,  by  replacing 
equations  (4.14)  with  appropriate  correlations  for  turbulent  heat  transfer  and 
skin  friction. 

If  the  thermodynamic  optimization  of  the  plate  fin  is  subjected  to  the 

constraint  of  fixed  material  (volume)  V  -  Re  Re  Re,  ,  the  entropy  generation 

L  b  o 

number  becomes 


N„  = 


.JUl 


1.15  V4  tanh(mh) 


+  1.328  B 


5  5 


(4.19) 


where  i  Re  /Re,  ,  n  =  Re  /Re  and  mL  =  1.15  y^^/(M  n .  Figure 

L  b  ob 

4.5  shows  the  trade-off  roll  played  by  the  geometric  aspect  ratio  y=  L/b  . 

The  optimum  aspect  ratio  y  for  minimum  irreversibility  can  be  determined 

opt 

numerically  based  on  calculations  of  the  type  illustrated  in  the  left  half  of 

Fig.  4.5.  The  general  behavior  of  the  y  ^  result  is  shown  in  the  right  half 

of  the  figure,  for  fixed  V  ,  M  arid  n.  The  optimum  aspect  ratio  y  reaches 

opt 

a  plateau  below  a  certain  value  of  friction  parameter  B  ;  as  B  increases 

_  ^ 

above  '  lu  ,  the  optimum  aspect  ratio  y  decreases  rapidly. 

opt 


Rectangular  plate  fin  of  trapezoidal  profile.  A  relatively  more  complex 
plate  fin  geometry  is  represented  in  Fig.  4.2(c).  This  time  the  fin  longitudinal 
section  is  trapezoidal,  again  in  the  spirit  of  E.  Schmidt's  principle  of  material 
(volume)  minimization  [4.5].  This  geometry  has  another  important  advantage  over 
the  rectangular  shape  of  Fig.  4.2(b):  the  tapered  profile  makes  trapezoidal 


optimization  (V  =  10  ,  M  =  100,  0  =  0.1).  Left  side:  entropy  generation  number  vs. 

:ct  ratio.  Right  side:  optimum  plate  aspect  ratio  vs.  friction  irreversibility  parameter 
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f  ins  accessible  to  metal-cutting  operat  ions  designed  to  remove  the  f  m-to-i  in 

material.  The  minimum  entropy  generation  design  of  this  class  of  fins  can  be 

approached  along  the  same  lines  as  the  design  of  rectangular  fins. 

Tn  the  interest  of  brevity,  we  illustrate  only  a  set  of  results  abstracted 

from  [4.l].  The  entropy  generation  rate  depends  on  four  geometric  parameters 

Re  ,  Re  ,  He,  ,  and  a  .  Whether  or  not  all  four  parameters  may  be  considered 
L  b  0  ‘ 

e 

independent  depends  on  specific  circumstances,  particularly  on  the  constraints 

faced  by  the  designer  (fixed  cost,  volume,  base  area,  etc.).  In  Fig.  4.6,  we 

show  numerical  results  obtained  in  the  case  when  b  ,  6  ,  and  a  are  fixed 

e 

by  design:  the  only  geometric  variable  in  this  case  is  the  fin  length  L 

(dimension  perpendicular  to  the  wall) .  Choosing  the  optimum  Re^  for  minimum 

entropy  generation  is  geometrically  equivalent  to  choosing  the  optimum  base 

thickness  Re..  since  Re.  =  Re.  +  2  Re  tan  a  .  The  left  side  of  Fig.  4.6 
0  0  0  L 

B  Be 

shows  that  N  has  a  sharp  minimum  with  respect  to  Re,.  .  The  optimum  values 
S 

of  base  thickness  Reynolds  number  are  reported  in  the  right  half  of  Fig.  4.6 

for  the  case  Re.  =  100  in  the  range  5°  <  a  <  10°  .  We  find  that  the 
o  —  — 

e 

optimum  fin  size  (Re^  )  decreases  as  the  fluid  friction  effect  (B)  becomes  more 

B 

pronounced,  in  agreement  with  conclusions  reached  in  previous  examples.  Further, 
wo  see  that  in  the  a  range  considered,  the  angle  a  has  a  relatively  minor 
impact  on  the  optimum  fin  size  for  minimum  irreversibility. 

Triangular  plate  fin  of  rectangular  profile.  To  the  thermodynamic  designer, 
this  geometry  is  challenging  due  to  the  absence  of  convenient  correlations  for 
heat  transfer  and  fluid  friction  in  the  three-dimensional  boundary  layer  flow 
which,  in  most  certainty,  will  cover  the  triangular  faces  of  the  fin.  This 
analytical  difficulty  can  be  partially  dealt  with  in  the  limit  D  <<  L  ,  where 
the  three-dimensional  effects  will  be  minor.  In  this  limit,  we  can  approximately 
treat  the  sharp-pointed  triangular  plate  as  a  flat  plate  in  parallel  flow,  with 


Arn 
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t he  special  property  that  the  length  swept  by  the  flow  (y)  is  a  function  of 
longitudinal  position  (x)  .  It  is  worth  pointing  out  that  the  D  L  limit 
is  in  full  agreement  with  the  uni-directional  heat  conduction  model  adopted 
in  sect  ion  4.1. 

The  analytical  path  leading  to  the  entropy  generation  rate  formula  is 
similar  to  the  method  used  in  the  earlier  examples  [4.l].  Again,  in  the 
interest  of  brevity,  we  report  only  a  number  of  representative  results. 

The  entropy  generation  rate  depends  on  three  geometric  parameters,  Re^  , 

Rej  and  Re,.  .  In  Fig.  4.7,  we  report  a  sample  of  optimum  fin  size  results 

for  cases  where  the  triangle  aspect  ratio  a  =  Re  /Re  is  fixed.  The  fin 

L  D 

irreversibility  reaches  a  clear  minimum  at  a  specific  value  of  fin  base 

width,  Re^  ;  the  optimum  fin  size  Re^  depends  on  the  relative  importance 
of  fluid  friction  irreversibility  (B) ,  on  the  metal-fluid  combination  (M) ,  and 
on  tiie  plate  thickness  (Re^).  The  right  hand  side  of  Fig.  4.7  is  a  summary 
ot  minimum  results  obtained  for  a  number  of  common  metal-fluid  combinations 

involving  copper,  aluminum,  water  and  air.  Regardless  of  combination,  the 


optimum  fin  size  (Re  ) 

D,opt 


decreases  as  the  triangular  shape  of  the  fin  becomes, 


by  design,  more  slender. 


Optimum  Fin  Matrix  for  Minimum  Entropy  Generation 


The  focus  of  the  preceding  sections  was  on  individual  fins.  Thus,  we 
developed  concrete  means  for  selecting  the  geometric  parameters  which  allow  the 
fin  to  perform  its  prescribed  heat  transfer  job  with  minimum  entropy  generation. 
An  equally  important  design  approach  consists  of  focusing  on  an  assembly  of  fins 
and  finding  out  which  particular  association  (matrix,  relative  positioning) 
allows  the  assembly  to  perform  its  prescribed  heat  transfer  job  with  minimum 
irreversibility.  We  illustrate  this  approach  in  this  section.  From  the  outset. 


generation  number  vs.  base  width.  Right  side:  optimum  base  width  vs.  geometric  aspec 
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we  note  that  the  task  of  determining  the  optimum  thermodynamic  association  of 
tins  is  made  difficult  by  the  lack  of  specific  heat  transfer  and  friction 
data  on  fin  matrices. 

The  optimum  association  of  pin  fins  can  be  determined  concretely  in  the 
"isothermal  fin"  limit:  in  this  limit  the  fin  matrix  is  well-approximated  by 
a  tube  bank  in  crossflow,  for  which  Kays  and  London  [4.6]  report  extensive  heat 
transfer  and  tluid  friction  data.  Consider  a  two-dimensional  array  of  staggered 
: ins.  The  overall  dimensions,  H  and  L,  and  the  total  heat  transferred 
by  the  fin  association,  Q,  are  considered  fixed  by  design.  The  total  rate 
of  entropy  generation  in  the  arrangement  is 
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Parameter  N  is  the  total  number  of  fins  populating  the  HxL  area.  For  heat 
transfer  and  drag  coefficients  we  have  [4.6] 

-2/3  -0.4  -0.18 

h  =  C.  P  U  c  Pr  Re,^  ,  C  =  C.  Re^  ,  (4.22) 

h  00  p  D  D  f  D 

where  parameters  C,  and  C  are  available  graphically  as  functions  of 
t  h 

transversal  and  longitudinal  pitch,  X  and  X  . 

t  t 

The  entropy  generation  number  Ng  emerges  as  a  function  of  the  matrix 
(  <  ,  ,  ;.  )  arid  the  diameter  of  each  fin  (Re^)  •  The  left  side  of  Fig.  4.8  shows 
the  dependence  of  N  on  Re^  in  one  situation  where  the  matrix  is  fixed, 
ib’. n  ly,  there  is  an  optimum  Re^  which  minimizes  the  total  entropy  generation 
rate;  corresponding  to  this  optimum  diameter,  there  is  an  optimum  number  of 


i 


i 


effect  of  friction  irreversibility  (B)  on  the  optimum 
fin  array.  Left  side:  the  optimum  pin  fin  diameter, 
optimum  number  of  fins. 
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fins  N  to  be  planted  on  the  surface  under  consideration.  The  right  side  of 

Fiq.  4.8  summarizes  the  results  of  the  minimization  procedure,  for  the 

case  where  only  the  longitudinal  pitch  is  held  constant.  The  optimum  fin 

diameter  increases  with  the  increasing  transversal  pitch  X  ,  while  the 

corresponding  number  of  fins  decreases. 

The  impact  of  varying  the  friction  parameter  B  is  illustrated  in  Fig.  4.9. 

If  the  matrix  (X^,  X  )  is  held  fixed,  the  optimum  fin  diameter  increases  as 

the  friction  effect  becomes  more  pronounced.  At  the  same  time,  the  total  number 

of  fins  (N)  corresponding  to  Re  ,  decreases. 

^  D ,opt 
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Nomenclature 


A  cross-sectional  area 

H  fluid  friction  irreversibility  parameter 

0  drag  coefficient 

L'  skin  friction  coefficient 

r 

!■'  drag  force 

h  enthalpy 

h  heat  transfer  coefficient 

k  thermal  conductivity  of  fin  material 

m  mass  flowrate 

N  total  number  of  isothermal  pin  fins  (tubes) 

entropy  generation  number 
P  pressure 

Pr  Prandtl  number 

o  total  heat  transferred  by  the  fin  association 

g"  local  heat  flux  per  unit  area 

q  heat  flux  through  the  fin  base 

s  entropy 

entropy  generation  rate 

■fen 

T  temperature 

u  velocity 

temperature  excess  at  the  fin  base,  T  -  T 
!S  B  “> 

'  thermal  conductivity  of  fluid 

u  dynamic  viscosity 

v  kinematic  viscosity 


dens i ty 


-76- 


5.  CONCLUSIONS 

This  final  report  summar ized  the  main  results  and  conclusions  of  a  two- 
year  study  of  the  thermodynamic  implications  of  heat  transfer  augmentation 
[1.6-1.12].  In  section  2  we  reviewed  a  number  of  central  topics  in  convective 
heat  transfer,  from  the  point  of  view  of  heat  transfer  as  an  entropy-producing 
mechanism.  We  saw  the  close  relationship  which  exists  between  the  irreversi¬ 
bility  due  to  heat  transfer  and  the  irreversibility  due  to  fluid  friction  in  a 
■liven  convective  heat  transfer  configuration.  In  section  2  we  developed  a 
number  of  analytical  tools  designed  to  facilitate  the  numerical  calculation  of 
thermodynamic  irreversibility  in  actual  applications. 

In  section  3  we  analyzed  quantitatively  the  thermodynamic  impact  of  using 
a  proposed  heat  transfer  augmentation  technique.  We  showed  that  a  technique's 
potential  for  reducing  the  rate  of  entropy  generation  depends  strongly  on  the 
<y,  ■orating  condit  ions  of  the  heat  exchanger  passage  in  question.  Specifically, 
the  augmentation  entropy  generation  number  depends  on  the  operating  Reynolds 
number  Re  and  on  the  irreversibility  distribution  ratio  <f  (the  ratio  of 
fluid  friction  irreversibility  divided  by  the  irreversibility  associated  with 
imperfect  thermal  contact).  When  the  operating  conditions  of  a  certain  component 
are  specified  (Re, if),  it  is  possible  to  assess  the  relative  thermodynamic  merit 
of  various  heat  transfer  augmentation  techniques  proposed  for  the  same  applica¬ 
tion.  This  comparison  is  done  on  the  basis  of  plots  (Re,(j>)  as  shown  in 

rigs.  3.2—3. 14;  each  of  these  plots  is  the  property  of  a  well-defined  heat 
transfer  augmentation  technique. 

finally,  in  section  4  we  considered  the  entropy  generation  associated  with 
the  use  of  extended  surfaces  (fins)  in  cross-flow.  Wc  showed  that  the  sizing 
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( t nermodynamic  selection)  of  an  individual  fin  always  involves:  a  trade-off 
between  heat  transfer  and  fluid  friction  i i reversibilities.  In  the  cases 
out  Line. 1  m  this  report,  the  thermodynamic  optimization  of  individual  fins 
could  be  carried  out  based  on  Uuily  simple  analysis.  In  general,  as  demon¬ 
s'  rated  in  a  more  detailed  report  [l. 12] ,  the  thermodynamic  selection  of 
tinned  surfaces  can  be  pursued  numerically.  The  main  conclusion  of  this 
last  phase  or  our  study  is  that  it  is  possible  and,  indeed,  advantageous  to 
coiceptual ize  finned  surfaces  which  perform  their  prescribed  heat  transfer 
duties  and,  at  the  same  time,  generate  the  least  amount  of  entropy  in  the 
engineering  system  in  which  they  function. 

C'voraJ 1,  our  study  illustrates  the  importance  of  the  Second  Law  of  Thermo¬ 
dynamics  m  the  development  of  ef f icienr.  power  systems.  Our  study  shows  how 
to  evaluate  the  change  in  thermodynamic  irreversibility  (entropy  generation) 
associated  with  a  proposed  desiqn  change,  and  how  to  make  the  right  exergy- 
auving  decisions.  Although  this  final  report  is  or1;,  a  partial  review  of  the 
wink  made  available  through  individual  publications  [1.6-1. 12] ,  it  is  hoped 
J  an  alone  serve  as  starting  point  for  future  Second  Law  analyses  of  power 


